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We review recent advances on tlie tlieory of spin qubits in nanostructures. We focus on 
four selected topics. First, we show how to form spin qubits in the new and promising 
material graphene. Afterwards, we discuss spin relaxation and decoherence in quantum dots. 
In particular, we demonstrate how charge fluctations in the surrounding environment cause 
spin decay via spin-orbit coupling. We then turn to a brief overview of how one can use 
electron-dipole spin resonance (EDSR) to perform single spin rotations in quantum dots 
using an oscillating electric field. The final topic we cover is the spin-spin coupling via spin- 
orbit interaction which is an alternative to the usual spin-spin coupling via the Heisenberg 
exchange interaction. 
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1. Introduction 

Spin qubits in quantum dot nanostructures (and the underlying physics) is a rapidly 
evolving research area. After the original proposal^^ based on the electron spin in few electron 
quantum dots, a number of related research branches have emerged. One of which deals 
with alternative ways to form spin qubits in solid state devices. Here, the aim is to find 
spin qubit realizations that either couple only weakly to the environment or that are easy to 
manipulate. A popular example with a growing interest in the spin qubit community is two- 
spin qubits (where the spin qubit corresponds to singlet and triplet states of two electrons). ^"^^ 
Further examples are many-spin cluster qubits composed of antiferromagnetically-coupled 
spin chains^' and spin qubits in magnetic molecules.^' ^'^^ Another active research field is 
devoted to coherence properties of spin qubits. Here, different aspects of spin relaxation and 
spin dephasing have been quantitatively analyzed for different dissipation channels. The most 
dominant ones are spin-orbit interaction, coupling the spin to lattice vibrations^ ^"^^^ and 
other charge fluctuations^^^ as well as the hyperfine interaction of the electron spin with the 
surrounding nuclear spins. ^^"'^^^ 

The success of spin qubits in nanostructures is substantially due to the major experimental 
breakthroughs that have been achieved in recent years (for recent review articles on spin 
qubits see Refs. 22-24). After pioneering experiments on few electron quantum dots,^^"'^^^ a 
first step towards the realization of quantum computing with the spin of electrons in quantum 
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dots has been made in single-shot measurements of the electron spin.^^^^^^ Subsequently, a 
coherent two-qubit gate (the \/SWAP gate) has been realized. '^^ Recently, coherent single 
spin rotations have been demonstrated via electron spin resonance techniques using pulsed 
magnetic fields. Thus, all single- and two-qubit operations required for universal quantum 
computing have been realized in spin qubits based on the original idea.^^ However, the time 
scales needed to operate single-qubit gates in spin qubits hosted in lateral quantum dots in 
GaAs/AlGaAs heterostructures^''"^^) are still quite long as compared to the decoherence time 
^2 ~ 1 — lO/US. (Note that this is not the case for two-qubit operations which can be performed 
as fast as 180 ps for the VSWAP gate.^^) The ratio of the operation time and the decoherence 
time should be of the order of 10'^ to be able to do fault-tolerant quantum computing. The 
decoherence time is currently limited by the hyperfine interaction of the electron spin with 
the surrounding nuclear spins. ^^"'^^^ Therefore, it is desirable to form spin qubits in other 
materials where spin relaxation and spin decoherence are less efficient than in GaAs/AlGaAs 
heterostructures. Two examples (which have already been realized) are few-electron quantum 
dots in carbon nanotubes^^"^^-* as well as semiconductor nanowires.^*^'^^) We discuss below 
another interesting example, namely spin qubits in graphene, where spin relaxation and spin 
decoherence mechanisms are expected to be weaker than in GaAs-based devices (for recent 
review articles on graphene see Refs. 42-44). 

The article is organized as follows: In Sec. 2, we explain in detail our recent proposal to 
form spin qubits in graphene quantum dots. Many of the aspects discussed in that section 
equally apply to spin qubits based on carbon nanotube quantum dots. In Sec. 3, different 
aspects of spin relaxation and decoherence are reviewed. Afterwards, in Sec. 4, recent ideas to 
use EDSR to form single spin rotations in quantum dots are discussed. In Sec. 5, we show how 
spin-spin coupling can be achieved via spin-orbit interaction. Finally, we conclude in Sec. 6. 

2. Spin qubits in graphene quantum dots 

It is generally believed that carbon-based materials such as nanotubes or graphene are 
excellent candidates to form spin qubits in quantum dots. This is because spin-orbit coupling 
is weak in carbon (due to its relatively low atomic weight), and because natural carbon 
consists predominantly of the zero-spin isotope ^^C, for which the hyperfine interaction is 
absent. In this section, we review how to form spin qubits in graphene.^^^ A crucial require- 
ment to achieve this goal is to find quantum dot states where the usual valley degeneracy is 
lifted. We show that this problem can be avoided in quantum dots with so-called armchair 
boundaries. We furthermore show that spin qubits in graphene can not only be coupled (via 
Heisenberg exchange) between nearest neighbor quantum dots but also over long distances. 
This remarkable feature is a direct consequence of the Klein paradox being a distinct property 
of the quasi-relativistic spectrum of graphene.^^^ 

Two fundamental problems need to be overcome before graphene can be used to form spin 
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Fig. 1. (Color online) A ribbon of graphcnc with semi-conducting armchair boundaries is schemati- 
cally shown. Two barrier gates (blue) define the rectangular size of the quantum dot (with width 
W and length L). A back gate (red) allows one to shift the energy levels in the dot. 



qubits and to operate one or two of them in the standard way.^'^^^ (i) It is difficult to create 
a tunable quantum dot in graphene because of the absence of a gap in the spectrum. ^^'^"^ (ii) 
Due to the valley degeneracy that exists in graphene, it is non-trivial to form two-qubit 
gates using Heisenberg exchange coupling for spins in neighboring dots. Several attempts 
have been made to solve the problem (1)^3-57) (-v^ithout having problem (ii) in mind). We have 
recently proposed a setup which solves both problems (i) and (ii) at once.^^^ In particular, 
we assume semiconducting armchair boundary conditions to exist on two opposite edges of 
the sample. It is known that in such a device the valley degeneracy is lifted, ^^'^''^ which is 
the essential prerequisite for the appearance of Heisenberg exchange coupling for spins in 
tunnel-coupled quantum dots, and thus for the use of graphene dots for spin qubits. 

We now discuss bound-state solutions in the appropriate setup, which are required for 
a localized qubit. We first concentrate on a single quantum dot which is assumed to be 
rectangular with width W and length L, see Fig. 1. The basic idea of forming the dot is to 
take a ribbon of graphene with semiconducting armchair boundary conditions in x-direction 
and to electrically confine particles in y-direction. 

The low energy properties of electrons (with energy e with respect to the Dirac point) in 
such a setup are described by the 4x4 Dirac equation 

cTxdx + crydy 

where the electric gate potential is assumed to vary stepwise, //(y) = /igate in the dot region 
(where < y < L), and ^(y) = /^barrier in the barrier region (where y < ov y > L). 
In Eq. (1), ax and dy are Pauli matrices (denoting the sublattices in graphene). The four 



^ + f^{y)^ = e^, (1) 
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component spinor envelope wave function * = {'i''^'',^'^'',-^'^' ,-*B ) 

varies on scales 

large compared to the lattice spacing. Here, A and B refer to the two sublattices in the 
two-dimensional honeycomb lattice of carbon atoms, whereas K and K' refer to the vectors 
in reciprocal space corresponding to the two valleys in the bandstructure of graphene. The 
appropriate semiconducting armchair boundary conditions for such a wave function can be 
written as (a = A, B)^^^ 

*f)U=v^ = e^'^/'^n=w. (2) 

These boundary conditions couple the two valleys and are, thus, the reason why the valley 
degeneracy is lifted. It is well known that the boundary condition (2) yields the following 
quantization conditions for the wave vector kx = Qn in x-direction^^' 

qn = {n± l/3)Tr/W, n G Z. (3) 

The level spacing of the modes (3) can be estimated as Ae Hvtt/SW, which gives Ae ~ 
30meV, where we used that v ~ 10^ m/s and assumed a quantum dot width of about W ~ 
30 nm. Note that Eq. (3) also determines the energy gap for excitations as E'gap = 2hvqo. 
Therefore, this gap is of the order of 60 meV, which is unusually small for semiconductors. 
This is a unique feature of graphene that will allow for long-distance coupling of spin qubits 
as will be discussed below. 

We now present in more detail the ground-state solutions, i.e. n = in Eq. (3). The 
corresponding ground-state energy e can be expressed relative to the potential barrier fj, = 
A^baxrier in the regions y < and y > L as e = /Xbarrier i ^^(^o + k^)^^"^- Here, the ± sign refers 
to a conduction band (-I-) and a valence band (— ) solution to Eq. (1). For bound states to 
exist and to decay at y — ^ ±oo, we require that Hvqo > |e — A'barrierl) which implies that the 
wave vector ky = k in y-direction, given by 

k = i\J ql - ((e - lJ-ha.vvier) / hvf , (4) 

is purely imaginary. In the dot region (0 < y < L) , the wave vector k in y-direction is replaced 
by k, satisfying e = //gate i hv(qQ + ^^)^/^. Again the it sign refers to conduction and valence 
band solutions. (In the following, we focus on conduction band solutions to the problem.) In 
the energy window 

/^bcirrier I, (5) 
the bound state energies are given by the solutions of the transcendental equation 



tan(^L) = ^ """-^ '""Z' ' r7r\2 - (6) 



hvk^y {hvqoY - (g - /Ubarricr)^ 
(e - /^barrier) (g " /^gate) " {fivqof ' 

We show a set of solutions to Eq. (6) for a dot with aspect ratio q^L = ttL/SW = 5 in Fig. 2. 
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Fig. 2. (Color online) Bound-state solutions of a dot with aspect ratio qqL = ttL/3W = 5. The 
diagonal lines indicate the region in which bound-state solutions do exist given by Eq. (5). All 
energies are taken in units of hvqQ. 

We now turn to the case of two coupled graphene quantum dots, separated by a potential 
barrier, each dot filled with a single electron. It is interesting to ask whether the spins Sj of 
these two electrons (i = 1,2) are coupled through an exchange coupling, H^xch = J^i ■ S2, 
in the same way as for regular semiconductor quantum dots,^^-* because this coupling is, in 
combination with single-spin rotations, sufficient to generate all quantum gates required for 
universal quantum computation.^-* The exchange coupling is based on the Pauli exclusion 
principle which allows for electron hopping between the dots in the spin singlet state (with 
opposite spins) of two electrons, but not in a spin triplet (with parallel spins), thus leading 
to a singlet-triplet splitting (exchange energy) J. 

However, a singlet-triplet splitting J ^ only occurs if the triplet state with two electrons 
on the same dot in the ground state is forbidden, i.e., in the case of a single non- degenerate 
orbital level. This is a non-trivial requirement in a graphene structure, as in bulk graphene, 
there is a two-fold orbital degeneracy of states around the points K and K' in the first Brillouin 
zone. This valley degeneracy is lifted in our case of a ribbon with semiconducting armchair 
edges, and the ground-state solutions determined by Eq. (6) are in fact non-degenerate. ^■^^ 
The magnitude of the exchange coupling within a Hund-Mulliken model is^^^ J = {—Uh + 
(Ufj + 16t'jf)^^'^)/2 + V , where t is the tunneling (hopping) matrix element between the left and 
right dot, U is the on-site Coulomb energy, and V is the direct exchange from the long-range 
(inter-dot) Coulomb interaction. The symbols tn and Uh indicate that these quantities are 
renormalized from the bare values t and U by the inter-dot Coulomb interaction. 

For t <^ U and neglecting the long-ranged Coulomb part, this simplifies to the Hubbard 
model result J = At^ /U where t is the tunneling (hopping) matrix element between the left 
and right dot and U is the on-site Coulomb energy. In the regime of weak tunneling, we 
can estimate t k, e J ^\{x,y)'^ fi{x,y)dx dy, where ^L,R{x,y) = ^{x,y ± (d -|- L)/2) are the 
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ground-state spinor wave functions of the left and right dots and e is the single-particle ground 
state energy. Note that the overlap integral vanishes if the states on the left and right dot 
belong to different transverse quantum numbers 7^ 

For the ground state mode, we have ul = tir = 0, and the hopping matrix element can 
be estimated for d > L as 

t K ieaod^WdzQ^k exp(-d|A:|), (7) 

where ao and Sq are wave function amplitudes (with dimension 1/Iength), see Ref. 48 for more 
details. As expected, the exchange coupling decreases exponentially with the barrier thickness, 
the exponent given by the "forbidden" momentum k in the barrier, defined in Eq. (4). 

The values of t, U, and J can be estimated as follows. The tunneling matrix element t is 
a fraction of £ ~ 30meV (for a width ofW^ 30 nm), we obtain that t ~ 0.5 . . . 2.5 meV. The 
value for U depends on screening which we can assume to be relatively weak in graphene,^^^ 
thus, we estimate, e.g., U ~ lOmeV, and obtain J ~ 0.1 . . . 1.5 meV. (Note that this rough 
estimate would correspond to very fast switching times ~ ^/J ~ 1 . . . lOps for the VSWAP 
operation.) 

For the situation with more than two dots in a line, it turns out that we can couple any 
two of them with the others being decoupled by detuning. In Fig. 3, we illustrate the situation 
of three dots in a line where the left and the right dot are strongly coupled and the center dot 
is decoupled by detuning. The tunnel coupling of dot 1 and dot 3 is then achieved via Klein 
tunneling through the valence band of the two central barriers and the valence band of the 
center dot. It is important for the long-distance coupling that the exchange coupling of qubit 
1 and qubit 3 is primarily achieved via the valence band and not via the qubit level of the 
center dot - leaving the qubit state of dot 2 unchanged. Using the standard transition matrix 
approach, we can compare the transition rate of coupling dot 1 and dot 3 via the continuum 
of states in the valence band of the center dot (which we call Fvb) with the transition rate 
via the detuned qubit level of the center dot (which we call Fqb). We obtain for the ratio^^^ 

rvB/TQE ~ (L/W) ln(4A/i5gap), (8) 

where A 6 eV is the band width of graphene. Therefore, by increasing the aspect ratio L/W, 
it is possible to increase the rate Fvb with respect to Fqb- For L/W = 2 and -Egap ~ 60meV, 
we find that Tyb/^QB ~ 12, meaning that the qubit level in dot 2 is barely used to couple dot 
1 and dot 3. This is a unique feature of graphene quantum dots due to the small and highly 
symmetric band gap. 

3. Spin relaxation and decoherence in quantum dots 

Phase coherence of spins in quantum dots (QDs) is of central importance for spin-based 
quantum computation in the solid state, however, the mechanisms of spin decoherence for 
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Fig. 3. (Color online) The energy bands of a triple quantum dot setup are shown in which dot 1 
and dot 3 are strongly coupled via cotunneling processes through the valence bands of barrier 2, 
barrier 3, and dot 2. The center dot 2 is decoupled by detuning. The energy levels are chosen such 
that A£2 ^ Aei. The triple dot example illustrates that in a line of quantum dots, it is possible 
to strongly couple any two of them and decouple the others by detuning. This is a unique feature 
of graphene and cannot be achieved in semiconductors such as GaAs that have a much larger gap 
(after Ref.^^)) . 



extended and localized electrons are rather different. Different mechanisms of spin relax- 
ation in QDs have been considered, such as spin-phonon coupling via spin-orbit (SO) in- 
teraction^^"^'^' or hyperfine interaction,^^) and direct hyperfine couphng.^^'^^^^^^ In this 
section, we consider the spin relaxation and decoherence in quantum dots due to the couphng 
to phonons and charge fluctuations in the surrounding environment.^^' We show how SO 
interaction couples the electron spin to these types of fluctuations by deriving an effective 
Hamiltonian for the spin subspace. Throughout this section, we consider only the leading 
contribution to SO interaction in two dimensional systems (commonly called linear-in-p SO 
interaction) 

HsO = (3{-PxCrx + PyCTy) + a{px(Jy - PyOx), (9) 

where a and (3 are Rashba and Dresselhaus coefficients, respectively. Both coefficients have 
been measured recently in GaAs/InGaAs quantum wells using optical detection schemes. ^"^^ 
(Note that in this and in the following sections the Pauli matrices denote the electron spin 
whereas in the previous section it was the sublattice index of graphene.) Moreover, we assume 
that the temperature is the smallest energy scale in the system and the Zeeman energy is less 
than the orbital quantization in the QD, Aj^T <C Ez <C fv^Q. 

Phonon contribution - Lattice vibrations perturb the confining potential U{r) of the dot 
and these fluctuations couple to the electron spin in the QD via the spin-orbit interaction. 
At low temperatures, the effective Hamiltonian for the electron spin is given 

^cff = ^gfiB[B + 6B{t)]-cT, (10) 
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B[T] 

Fig. 4. (Color online) Solid curve: The relaxation rate 1/Ti due to phonons as a function of an in-plane 
B for a GaAs QD with huo = 1.1 meV, Xgo = fi/m*(3 = 1 fim, and a — 0. Dashed (dotted) curve: 
Contribution of the piezoelectric mechanism with transverse (longitudinal) phonons. Dot-dashed 
curve: Contribution of the deformation potential mechanism. Different coordinate frames are used 
for the relaxation rate calculations; {x,y,z) are the main crystallographic axes and {x',y',z') are 
defined as x' — {x + y)/V2, v' — {u — and z' = z. 



6B{t) = 25 X n{t), (11) 

where B is the applied magnetic field and Q,{t) is the quantum fluctuating field due to the 
coupling to phonons. Eqs. (10,11) show an important result: In first order in SO interaction, 
there can be only transverse fluctuations of the effective magnetic field, i.e., 6B{t) ■ B = 0, 
and the coupling is proportional to the i?-field itself. The former property holds true for 
spin coupling to any fluctuations, be it the noise of a gate voltage or coupling to particle- 
hole excitations in a Fermi sea. Consequently, there is no pure dephasing and T2 = 2Ti for 
arbitrarily large Zeeman splitting, in contrast to the naively expected case T2 <^ Ti, where 
Ti is the longitudinal relaxation time (or simply relaxation time) and T2 is the transverse 
relaxation time (decoherence time) of the spin. After averaging over the phonon bath, we 
find that the spin decay rate has a non-trivial magnetic field dependence; we do not present 
here the full analytic expression for the magnetic field dependence of the spin decay rate and 
refer to the original work instead. ^'^^ We plot 1/Ti as a function of in-plane B for a = 
(only Dresselhaus spin-orbit interaction), see Fig. (4). In agreement with experiment, 
shows a plateau in a wide range of B fields, due to a crossover from piezoelectric-transverse 
(dashed curve) to the deformation potential (dot-dashed curve) mechanism of electron-phonon 
interaction. Note that if a = /3 and B \\ y' then 1/Ti vanishes (the same is true for a = —f3 and 
B II x'), where x' = [110] and y' = [110]. A detailed measurement of the -B-field dependence 
was reported recently giving very good agreement with theory. Quite remarkably, the 
largest measured Ti times exceed 1 second. 
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Quantum point contact (QPC) contribution - Charge fluctuations in the surrounding en- 
vironment of the QD cause spin decay. Here, we consider one of these sources, a nearby 
functioning QPC (see Fig. 5), in which the charge couples to the spin via spin-orbit inter- 
action in the presence of a magnetic field. The effective Hamiltonian for the electron spin 
looks the same as Eq.(lO), but here, the origin of Vt[t) is the electron shot noise in the QPC 
and its functional dependence on the system parameters is different from the phonon case. 
There are two mechanisms which contribute to the spin relaxation rate 1/Ti: The electron- 
hole excitations in the QPC Fermi leads and the electron shot noise in the QPC.^^^ In the 
regime with high bias voltages AyU applied to the QPC, the latter is the dominant one. To go 
further, we assume that the applied magnetic field B is in-plane and along x' and we obtain 
{Ez, r < |A// ± Ez\ < hivoY^^ 

Sll = ^r(i-r). (13) 

Here u = 1/2'iThvF is the density of states per spin and mode in the QPC leads, m* is the 
electron effective mass, n is the dielectric constant, a is the distance form the QD center to the 
QPC, is the orientation angle of the QPC on the substrate (see Fig.5), A^c is the Coulomb 
screening length, A± = h/m*{(3 it a) are spin-orbit lengths, Jlwo is the orbital quantization 
energy in the QD, T is the transmission coefficient of the QPC, and Sll is the current shot 
noise. Therefore, in this regime, spin decay rate is linear in bias voltage A/x and scales as a"^. 
Moreover, Ti strongly depends on the QPC orientation on the substrate (the angle 9 between 
the axes x' and X, see Fig. (5)), e.g. the non-equilibrium part of the relaxation rate vanishes 
at = 7r/2, for an in-plane magnetic field B along x'. We conclude that the spin decay rate 
can be minimized by tuning certain geometrical parameters of the setup. Our results should 
also be useful for designing experimental setups such that the spin decoherence can be made 
negligibly small while charge detection with the QPC is still efficient. 

4. EDSR in quantum dots 

Spin-orbit interaction, although it is one of the main sources of the spin decay in QDs, 
can be employed to manipulate the electron spin. Here we show how by using an ac electric 
field together with a static magnetic field, one can coherently rotate the spin of the electron 
around the Bloch sphere.^^^ The physical mechanism responsible for the spin rotation is the 
so-called EDSR^^^ which is similar to usual Electron Spin Resonance (ESR)^^'^^) but, in the 
former case, the oscillating electric field replaces the oscillating magnetic field in the latter 
case. The main advantage of EDSR to ESR is its experimental convenience. 

We derive an effective Hamiltonian for the electron spin in the presence of a coherent 
driving ac electric field V{r,t) = e dr' ■ E{r',t) eE{t) ■ r (see Fig. 6). Hereby, we use 
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Fig. 5. Schematic of the quantum dot (QD) coupled to a QPC. The {X,Y) frame gives the setup 
orientation, left (L) and right (R) leads, with respect to the crystallographic directions x' = [110] 
and y' = [110]. The dot has a radius Ad and is located at a distance a from the QPC. The vector 
R describes the QPC electrons and r refers to the coordinate of the electron in the dot. The noise 
of the QPC current / perturbs the electron spin on the dot via the spin-orbit interaction. 



the dipole approximation, ignoring the coordinate dependence of the electric field due to 
the smallness of the dot size compared to the electric field wavelength. This leads us to the 
following effective spin Hamiltonian^^^ 

^eff = ^gfiBB-<T+^h{t)-(T, (14) 

h{t) = 2gfiBBxn{t), (15) 
= -^{XZ'Ey,{t),\^'E,,{t),0) , (16) 

where E{t) = Eq sin {uact) of amplitude Eq = Eo{cos sin cf), 0) and is the angle of Eq with 
respect to the axis x' (see Fig.6). Note that the resonance happens when u>ac = u>z = Ez/h, i.e. 
when the frequency of the driving field matches the Larmor frequency. The above Hamiltonian 
has a similar form to the ESR Hamiltonian, except for the fact that the oscillating electric 
field plays the role of the ac magnetic field. Consequently, we can rotate the electron spin 
around the Bloch sphere and build a universal single qubit gate. However, to quantify the 
efficiency of our EDSR scheme, we need to estimate the amplitude of the EDSR field, h{t), 
which is proportional to the Rabi frequency ujr. For GaAs QDs, we assume that A+ A_ 
Xso = 8 /xm, l^l = 0.44, hujQ = 1 meV, and Eo = 10^ V/cm, which yields ~ 10~^. Together 
with the applied magnetic field 5 = lOT, we obtain lor ~ 10^ s~^. We conclude that, with 
the present QD setups, EDSR enables one to manipulate the electron spin on a time scale of 
10 ns, which is considerably shorter than typical spin dephasing times T2 ~ 1 — lO/US in gated 
GaAs QDs. This mechanism has recently been employed to experimentally rotate the electron 
spin in quantum dots^^^ with 7r/2 rotations as fast as ~ 55ns. In a similar experimental setup, 
hyperfine-mediated gate-driven electron spin resonance has been observed. 
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Fig. 6. Setup for electric field control of spin via the spin-orbit interaction. The quantum dot (QD) 
contains a single electron with spin S = {h/2)a, deep in the Coulomb blockade valley, and in the 
presence of an external static magnetic field giving rise to a Zeeman splitting Ez- The gates 1 and 
2 are used to generate an alternating electric field E{t), which acts via the spin-orbit interaction 
on the electron spin. As a result, an electric dipole spin resonance (EDSR) occurs if the frequency 
of E{t) is tuned to match the Larmor frequency Wz = Ez/h. 



Up to now, we have only considered the linear-in-jo spin-orbit interaction. However, if the 
two dimensional electron gas (2DEG) has a finite width d, then the so-called terms of the 
Dresselhaus spin-orbit interaction also come into the play: 

HsO = 2 {PyPxPy(^x - PxPyPxf^y) , (17) 

where 7 = ad \/l'm*'^Eg is the spin-orbit coupling constant, with (~ 0.07 for GaAs), and 
Eg the band gap. Quite remarkably, if the quantum dot potential is harmonic, then the spin 
does not couple to E{t) in the first order of Hgo and zeroth order of Ez^^^ Thus, for a 
harmonic confining potential, one is left with the same dominant mechanism as considered 
above for the "linear in terms. To estimate the strength of the resulting EDSR, we expand 
in terms of the Zeeman interaction and note that 7 ~ (5(P If?^ and therefore the amplitude 
of h{t) = 2giiBB X Q,{t) is by a factor d^/A^ ^ 1 smaller as compared to the corresponding 
amplitude of the linear-in-p contributions. 

Next we consider a quantum dot with anharmonic potential U{r) and show that the p^- 
terms in Eq. (17) give rise to a spin-electric coupling proportional to the cyclotron frequency 
u)c = eBz/m*c.^^^ Since hujc differs parametrically from Ez {Ez/hwc = gm*B/2mBz), the 
p^-terms can be as significant as the p-terms, provided Ez/hujc ^ ^ l^i which is realistic for 
GaAs quantum dots. As an example, we consider i7(r) = frfuj^r^ jl -|- Tyr^ , where 77 is a 
measure of deformation from a harmonic confinement, and obtain^^^ 

1 G^Tlfi (jj 

2^^^^ ■ = {Ey{t)a, + E,{t)ay) . (18) 

Finally, we note that the p^-terms can also be relevant for spin relaxation in quantum dots 
with anharmonic confining potential. Of course, the magnetic field has to have an out-of-plane 
component for this spin-electric coupling to dominate over the one considered in the previous 
section. 
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In this part of the article, we discuss the interaction of two electron spins localized in 
quantum dots through the combined effect of spin-orbit interaction and Coulomb repulsion. 
The two single-electron quantum dot system is shown in Fig. 7. It is assumed that the two 
dots are well separated from each other such that there is no electron tunneling between 
them. In this respect, the interaction between the spins is fundamentally different from the 
Heisenberg exchange interaction for which the presence of tunneling is crucial. Similarly, 
the combined effect of Heisenberg exchange interaction and spin-orbit coupling^^"^^)is also 
based on tunneling and should be carefully distinguished from the spin-orbit effect studied 
here. Even though the Heisenberg exchange coupling allows typically for much stronger spin- 
spin coupling than the electrostatically induced one,^^) the latter one can prove useful for 
cases where it is difficult to get sufficient wavefunction overlap (needed for large Heisenberg 
exchange), and, moreover, it is also important to understand in detail the electrostatically 
induced spin-spin coupling in order to get control over possible interference effects between 
different types of coupling. This will be of importance for spin-qubit applications in order to 
minimize spin decoherence and gate errors. 

We give now a short theoretical description of our system. The Hamiltonian of the two- 
single electron quantum dot system is 

^ = E + ^(''^) + Idf^B^ ■ + H^so) + n — ^ P (19) 

.f^2 V^"^ ^ / '^1^1 - r2 -I- ao| 

where the first two terms are the kinetic and orbital confinement {U{ri) = m*u)Qrf /2), the 
third term is the Zeeman energy, the fourth term stands for the spin-orbit interaction, both 
Rashba and Dresselhaus [see Eq. (9)], while the last term stands for the Coulomb coupling 
between the two electrons. The distance between the centers of the two dots is oq. Usually, 
the spin-orbit interaction is a weak perturbation compared with the orbital level spacing and 
as a consequence can be treated within perturbation theory, as it was done also in the two 
previous sections. However, here we have an additional energy scale given by the strength of 
the Coulomb repulsion. This strength is measured through the parameter 5 = {X/aB)iX/ao)J^^ 
where A = ^JhfmFujQ is the dot radius and an = h^n/m*e^ is the Bohr radius in the material. 
In the general case of arbitrary strong Coulomb repulsion, the effective spin Hamiltonian Hspin 
of the two electron system reads^^^ 

TT 1 ij^eif 1 I 1 ij^eif 2 I T 12,7 12 /on\ 

iispin = 2^iz<^z + 2^2Z^z + ^xf^x^x + Jy^^y^^y, l^Uj 

where the explicit expressions for the spin-orbit renormalized Zeeman splittings E'^^ and 
spin-spin couplings Jx,y are given in Ref. 76. This interaction vanishes for vanishing Zeeman 
splitting and is highly anisotropic (XY type). Due to the finite Zeeman splitting, the relevant 
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Fig. 7. (Color online) The figure shows a sketch of the model system which consists of two identical 
quantum dots in the xy-plane, separated by distance ao (measured from dot-center to dot-center). 
Si denotes the spin of electron i = 1, 2, A is the dot radius, and B is the external magnetic field. The 
respective orbital wave functions of electron 1 and 2 are assumed to have no overlap (i.e. tunneling 
between the dots is excluded). The remaining purely electrostatic Coulomb interaction between 
the electron charges leads, via spin-orbit interaction, to an effective coupling between their spins. 
This spin-spin interaction depends sensitively on the orientation of B, with no component along 
it, and is proportional to B^. 



electrostatically induced spin-spin coupling can be written as 

Hs-s = Jeff{<yW- + (21) 

with cr± = (Tx i iay and Jeff = (1/2)(J^ -|- Jy)-'^^^ Up to now we posed no assumptions on the 
strength of the Coulomb repulsion. However, there are two interesting limiting cases, namely 
(5 <C 1 (weak Coulomb repulsion) and (strong Coulomb repulsion). 

In the first case, 5 <C 1, the Coulomb interaction is a weak perturbation compared to the 
bare orbital level spacing hwQ such that 

Hs-s = I drrdr, , ^^^^^^ , . (22) 
J K\ri-r2 + ao\ 

Here, the 2x2 matrices 5pi^2 are the spin-orbit induced charge distributions or spin-dependent 

charge distributions in each dot in the absence of Coulomb interaction.^^) From Eq. (22) we 

see that the spin-spin interaction results from a Coulomb-type coupling between two charge 

distributions which themselves depend on spin. In the limit of large interdot distances oq ^ A, 

we can perform a multipolar expansion, such that within the lowest order we obtain 

mi ■ m2 - ^{mi ■ na){m2 ■ ria) 

Hs-s ~ 5 (^•Jj 

where = aQ/uQ. The dipole moments rrij = {0\6piri\0) = fiai, where |0) is the orbital 
ground-state and /2 is the tensor corresponding to an effective spin-orbit magneton (for explicit 
expressions see Ref . 76) . The strength of this effective spin-orbit induced magneton is given by 
||/i|| PS eEz /m*ujQ\so- To give an estimate, we assume fuujQ ^ O.SmeV, Ez ~ 0.05 meV {B ~ 
2T) and m* = 0.067me, Xso ~ 10^® m for GaAs quantum dots which gives, when compared 
with the Bohr magneton, ~ 10'^. This implies that the spin-orbit induced dipole- 

dipole interaction in Eq. (23) can be much stronger than the direct dipole-dipole interaction 
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Fig. 8. (Color online) The function G occurring in Eq. (24) plotted as a function of the geometric 
distance qq between the dot centers scaled by the dot radius A for different magnetic field orien- 
tations. The dashed line represents the dipolar approximation of G for a perpendicular magnetic 
field (6 — 0) which scales like Gq^. 



in vacuum, whose strength is given hy fiB- Also, still in the limit (5 <C 1, but for arbitrary 
interdot distance oq the effective coupling Jf,ff has the form 

where the function G(ao/A, 9, ^) is plotted in Fig. 8 as a function of a^/X for different angles 
9,^. The key feature of the electrostatic spin-spin interaction is that it can range from fer- 
romagnetic to antiferromagnetic type, depending on the magnetic field orientation, passing 
even through zero for certain angles and/or inter-dot distances. 

We now focus on the opposite limit S ^ 1, when the Coulomb interaction is much stronger 
than the bare orbital level spacing Hluq. Then, we approximate e^/K|r-|-ao| — > (e^/2Ka^)[3(na • 
r)^ — r^],''^) where a is the effective distance between the electrons due to the combined effect 
of Coulomb repulsion and orbital confinement. For the explicit derivation of the effective 
distance a in terms of the bare one qq see Ref. 76. Within this ansatz, the spin-spin coupling 
Hamiltonian takes the form 



El 



(25) 



for the case of a perpendicular magnetic field. In the above Gxprcssion wg liavG bx — 
Vl + 4(A/os)(A/a)3 and by = y^l - 2(A/aB)(A/a)3. 

Let us give now some estimates for the coupling J^jj when an in-plane magnetic field is 
applied along, say, the x-direction. Assuming now GaAs quantum dots, and Ez = 0.1 meV 
{B = 4 T), huJo = 0.5 meV (A/a^ ~ 5), X/Xso ~ 10^^. Using these numbers and taking for 
the geometric inter-dot distance ao/A ~ 2, we obtain J^jj ~ 10"'' eV. It is worth mentioning 
that the hyperfine interaction between the electron and the collection of nuclei in a quantum 
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dot (~ 10^) leads to similar energy scales. ^^'^'^^ This shows that the spin-spin coupling derived 
here can be very relevant for the spin dynamics in the case of electrostatically coupled quantum 
dots and that it can also compete with other types of interactions. Considering now the case 
of InAs quantum dots'^'^'^^^ in a magnetic field along the x direction, with A50 ~ 2A ~ lOOnm 
and Ez/Hloo = 0.1 and taking also oq/A ~ 2, a value of J^ff ~ 10~^eV is obtained. 

6. Conclusions 

We have discussed several selected topics on the theory of spin qubits in nanostructures. 
We have first reviewed our recent proposal how to form spin qubits in graphene. This is 
interesting for two reasons. On the one hand, one expects very long spin lifetimes in graphene 
because of a weak spin-orbit interaction and very few host atoms with a nuclear spin. On the 
other hand, spin qubits in graphene allow for a new type of long distance coupling that uses the 
property that a ribbon of graphene is a small bandgap semiconductor. Furthermore, we have 
pointed out several aspects of spin relaxation and decoherence due to spin-orbit interaction 
and the coupling to a bath. As two possible dissipation channels we have considered lattice 
vibrations (phonons) and charge fluctuations in the surrounding environment, for instance, a 
nearby quantum point contact. Subsequently, we have shown how to use EDSR to rotate the 
spin of an electron in a quantum dot using an oscillating electric field (instead of the oscillating 
magnetic field employed in the usual ESR). In the final part of the review article, we have 
discussed how to couple two spins (located in two different quantum dots) via spin-orbit 
interaction in a situation in which direct tunneling between the dots is highly suppressed. 

We would like to thank D.V. Bulaev, G. Burkard, and V.N. Golovach for the collaboration 
on the work reviewed in this article. Financial support has been provided by the Swiss NSF, 
the NCCR Nanoscience, and JST ICORP. 
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